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Abstract. This paper is concerned with the distributed optimal control of a time-discrete Cahn- 
Hilliard/Navier-Stokes system with variable densities. It focuses on the double-obstacle potential 
which yields an optimal control problem for a family of coupled systems in each time instance of 
a variational inequality of fourth order and the Navier-Stokes equation. By proposing a suitable 
time-discretization, energy estimates are proved and the existence of solutions to the primal system 
and of optimal controls is established for the original problem as well as for a family of regularized 
problems. The latter correspond to Moreau-Yosida type approximations of the double-obstacle 
potential. The consistency of these approximations is shown and first order optimality conditions for 
the regularized problems are derived. Through a limit process, a stationarity system for the original 
problem is established which is related to a function space version of C-stationarity. 
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1. Introduction. In this paper we are concerned with the optimal control of 
two (or more) immiscible fluids with non-matched densities. For the mathematical 
formulation of the fluid phases, we use phase field models which have recently been 
used successfully in applications involving, e.g., phase separation phenomena (see, e.g., 
[4, 14, 30]). Some of the strengths of phase field approaches are due to their ability 
to overcome both, analytical difficulties of topological changes, such as, e.g., droplet 
break-ups or the coalescence of interfaces, and numerical challenges in capturing the 
interface dynamics. In this context, a so-called order parameter depicts the concentra¬ 
tion of the fluids, attaining extreme values at the pure phases and intermediate values 
within a thin (diffuse) interface layer, and it is associated with decreasing/minimizing 
a suitably chosen energy. 

A renowned diffuse interface model is the Cahn-Hilliard system which was first 
introduced by Cahn and Hilliard in [9]. In the presence of hydrodynamic effects, the 
system has to be enhanced by an equation which captures the behavior of the fluid. 
In [28], Hohenberg and Halperin published a first basic model for immiscible, viscous 
two-phase flows. Their so-called ’model H’ combines the Cahn-Hilliard system with 
the Navier-Stokes equation. It is however restricted to the case where the two fluids 
possess nearly identical densities, i.e., matched densities. Recently, Abels, Garcke 
and Grim [2] obtained the following diffuse interface model for two-phase flows with 
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non-matched densities: 


(l.la) 

(1.1b) 

(1.1c) 

(l.ld) 

(l.le) 

(l.lf) 

(1-lg) 


dt<p + v\7tp — div(m(v?)V/i) 
—A tp + d^o (tp) — p — ntp 
d t (p(tp)v) + div(w <g> p(tp)v) — div(2p(tp)e(v)) + X7p 

+div(w <g)«/) — pS7tp 
divv 


v\on 

9 n p\dn = d n p\an 

(v,<p)\t=o 


0 , 

0, 

0, 

0, 

0, 

0 , 

(Va,<Pa), 


which is supposed to hold in the space-time cylinder fl x (0,oo), where dfl denotes 
the boundary of f 1. This system is thermodynamically consistent in the sense that it 
allows for the derivation of local entropy or free energy inequalities. 

In the above model, v represents the velocity of the fluid and p describes the 
fluid pressure. The symmetric gradient of v is defined by e(v) := f (Vu + Vu T ). The 
density p of the mixture of the fluids depends on the order parameter tp which reflects 
the mass concentration of the fluid phases. More precisely, 

(1-2) PVP) = -j- + -2-^ 


where tp ranges in the interval [—1,1], and 0 < p\ < p 2 are the given densities of 
the two fluids under consideration. This is one of the main distinctions compared 
to the model with matched densities, where p is a fixed constant. The quantity p 
denotes the chemical potential in the Cahn-Hilliard system and helps to split the 
fourth-order in space differential operator into two second-order operators. Another 
important difference between (1.1) and model ’H’ is the presence of a relative flux J := 
— P 2 "? 1 m(tp)Vp which corresponds to the diffusion of the two phases and additionally 
complicates the analytical situation. The viscosity and mobility coefficients of the 
system, rj and m , depend on the actual concentration of the two fluids at each point 
in time and space. The initial states are given by v a and tp a , and k > 0 is a positive 
constant. Furthermore, 'Fo represents the convex part of the homogeneous free energy 
density contained in the Ginzburg-Landau energy model which is associated with the 
Cahn-Hilliard part of (1.1). Usually, the homogeneous free energy density serves 
the purpose of restricting the order parameter tp to the physically meaningful range 
[—1,1] and to capture the spinodal decomposition of the phases. For this reason, it 
is typically non-convex and maintains two local minima near or at —1 and 1. 

Depending on the underlying applications, different choices have been investi¬ 
gated in the literature. In their original paper [9], Cahn and Hilliard considered the 
logarithmic form 'F(vj) = (1 + y>)ln(l + tp) + (1 — tp) ln(l — tp) — ^ tp 2 which also 
plays an important role in the Flory-Huggins solution theory of the thermodynamics 
of polymer solutions. Another possible choice is the smooth double-well potential 
^{tp) = f (1 — tp 2 ) 2 , see e.g. [12, 16]. It permits pure phases but fails to restrict 
the order parameter to [—1,1]. Therefore, it is perhaps a less relevant choice in ma¬ 
terial science. In [35], Oono and Puri found that in the case of deep quenches of, 
e.g. binary alloys, the double-obstacle potential, is better suited than the other free 
energy models mentioned above. A similar observation appears to be true in the case 
of polymeric membrane formation under rapid wall hardening. The double-obstacle 
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potential ^(tp) = Fr-iu^) — with denoting the indicator function of 

the interval [—1,1] in R, combines the advantage of the existence of pure phases and 
the exclusiveness of the interval [—1,1] at the cost of losing differentiability (when 
compared, e.g., to the double-well potential). As a consequence, (1.1b) becomes a 
variational inequality which complicates the analytical and numerical treatment of 
the overall model. 

In this paper we study the optimal control of a time discrete coupled Cahn- 
Hilliard-Navier-Stokes (CHNS) system with the double-obstacle potential. For this 
purpose, we introduce a distributed control u which enters the Navier-Stokes equation 
(1.1c) on the right-hand side and aims to minimize an objective functional J subject 
to the control-version CHNS(u) of the Cahn-Hilliard-Navier-Stokes system: 

minimize J {p, p,v,u) over (p, p, v, u ) 

subject to (s.t.) u € Uad, (</>,satisfies CHNS(u), 

where U a d is a given set of admissible controls. 

Regarding physical applications, we point out that the CHNS system is used to 
model a variety of situations. These range from the aforementioned solidification 
process of liquid metal alloys, cf. [14], or the simulation of bubble dynamics, as in 
Taylor flows [4], or pinch-offs of liquid-liquid jets [29], to the formation of polymeric 
membranes [45] or proteins crystallization, see e.g. [30] and references within. Fur¬ 
thermore, the model can be easily adapted to include the effects of surfactants such 
as colloid particles at fluid-fluid interfaces in gels and emulsions used in food, phar¬ 
maceutical, cosmetic, or petroleum industries [5, 37]. In many of these situations an 
optimal control context is desirable where the system is influenced in such a way that 
a prescribed system behavior needs to be guaranteed. 

In the literature, the classical case of two-phase flows of liquids with matched 
densities is well investigated, see e.g. [28]. When it comes to the modeling of fluids 
with different densities, then the literature presents various approaches, ranging from 
quasi-incompressible models with non-divergence free velocity fields, see e.g. [32], to 
possibly thermodynamically inconsistent models with solenoidal fluid velocities, cf. 
[13]. We refer to [7, 8, 18] for additional analytical and numerical results for some of 
these models. In [1], Abels, Depner and Garcke derived an existence result for the 
given system (1.1) with a logarithmic potential, and in the recent preprint [19] system 
(1.1) with smooth potentials (thus excluding the double-obstacle homogeneous free 
energy density) is considered in a fully discrete and an alternative semi-discrete in 
time setting including numerical simulations. 

The optimal control problem associated to the Cahn-Hilliard-Navier-Stokes sys¬ 
tem with matched densities and a non-smooth homogeneous free energy density 
(double-obstacle potential) has been previously studied by the first and last author of 
this work in [27]. We also mention the recent preprint [17] which treats the control of 
a nonlocal Cahn-Hilliard-Navier-Stokes system in two dimensions. Apart from these 
contributions the literature on the optimal control of the coupled CHNS-system with 
non-matched densities is - to the best of our knowledge - essentially void. Neverthe¬ 
less, we mention that there are numerous publications concerning the optimal control 
of the phase separation process itself, i.e. the distinct Cahn-Hilliard system, see e.g. 
[10, 11, 21, 25, 43, 44], 

We point out that the presence of a non-smooth homogeneous free energy density 
associated with the underlying Ginzburg-Landau energy in the Cahn-Hilliard system 
gives rise to an optimal control problem for the Navier-Stokes system coupled to 
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the Cahn-Hilliard variational inequality. In particular, due to the presence of the 
variational inequality constraint, classical constraint qualifications (see, e.g., [46]) fail 
which prevents the application of Karush-Kuhn-Tucker (KKT) theory in Banach space 
for the first-order characterization of an optimal solution by (Lagrange) multipliers. In 
fact, it is known [22, 27] that the resulting problem falls into the realm of mathematical 
programs with equilibrium constraints (MPECs) in function space. Even in finite 
dimensions, this problem class is well-known for its constraint degeneracy [33, 36]. 
As a result, stationarity conditions for this problem class are no longer unique (in 
contrast to KKT conditions); compare [22, 23] in function space and, e.g., [39] in 
finite dimensions. They rather depend on the underlying problem structure and/or 
on the chosen analytical approach. In this work, we utilize a Yosida regularization 
technique with a subsequent passage to the limit with the Yosida parameter in order 
to derive conditions of C-stationarity type. This technique is reminiscent of the one 
pioneered by Barbu in [6], but for different problem classes. 

The remainder of the paper is organized as follows. In section 2 we introduce the 
semi-discrete Cahn-Hilliard-Navier-Stokes system and assign it to the corresponding 
optimal control problem. In section 3 we show the existence of feasible points to 
the original optimal control problem, as well as to regularized problems. Section 4 is 
concerned with the existence of globally optimal solutions, and section 5 deals with 
the consistency of the chosen regularization technique. In section 6 we derive first- 
order optimality conditions for the regularized problems using a classical result from 
non-linear optimization theory. Then a limiting process leads to a stationarity system 
for the original problem. The latter is the content of section 7. 

2. The semi-discrete CHNS-system and the optimal control problem. 

As a first step towards the numerical treatment of the underlying Cahn-Hilliard- 
Navier-Stokes system, we study a semi-discrete (in time) variant. For our subsequent 
analysis we start by fixing the associated function spaces and by invoking our working 
assumptions. 

For this purpose, let C R n ,N = 2,3, be a bounded domain with smooth 
boundary dfl £ C 2 (f2). In particular, f2 satisfies the cone condition, cf. [3, Chapter 
IV, 4.3], 

For k € N and 1 < p < oo we introduce the following Sobolev spaces: 

Hq R n ) = {/ g H k (Cl-R N )nH^(Cl-R N ) : div/ = 0, a.e. on fi} , 

W k ’ P (n) = | f £ W k ’ p (n) : J fdx = 0 

W^(O) = {/ G W fc ’ P (H) : d n f [d n = 0 on 512 } , 

where ’a.e.’ stands for ’almost everywhere’. Here, W k ’ p (fl) and W"Q’ p (f2) denote 
the usual Sobolev space, see [3]. For p = 2, we also write H k (Q) respectively Hq (f2) 
instead. Unless otherwise noted, (•, •) represents the L 2 -inner product, ||-|| the induced 
norm, and (•,•} := (•,-)-^-i -^-i the duality pairing between H (f2) and H (f2). For 
a Banach space W, we denote by W* its topological dual, and C(W, W *) defines the 
space of all linear and continuous operators from W to W* . In our notation for norms, 
we do not distinguish between scalar- or vector-valued functions. The inner product 
of vectors is denoted by ’•’, the vector product is represented be ’(g)’ and the tensor 
product for matrices is written as 
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Remark 2.1. Before we present the semi-discrete system and assuming integra- 
bility in time, from (1.1a) we get 

/ dtpdx = — / v\7pdx + / div(m(<p)Vfi)dx = 0, 

Jq J q, Jq 

Hence utilizing (1-lg) the integral mean of p satisfies 



i.e., it is constant in time. By assuming p a £ (—1,1), we exclude the uninteresting 

case \Tp^\ = 1. This can be achieved by considering the shifted system (1.1), where p 

_2 

is replaced by its projection onto L (f2). Consequently, we need to work with shifted 
variables such as, e.g. m(y + Tpf), which we again denote by m(y) in a slight misuse 
of notation. 

Motivated by physics, we assume throughout that the mobility and viscosity 
coefficients are strictly positive as specified in Assumption 2.2 below. Furthermore, 
we extend the connection (1.2) between p and p to all of R, as our studies include 
certain double-well type potentials which allow for values of p outside the physically 
relevant interval [—1,1]. 

Assumption 2.2. 

1. The coefficient functions m,r) £ C 2 (R) in (1.1c) and (1.1a) as well as their 
derivatives up to second order are bounded, i. e. there exist constants 0 < bi < 
b 2 such that for every x £ R, it holds that b\ < min{m(x), 7/(x)} and 

max{m(x), r/(x), |m'(x)|, \r}'(x)\, |m"(x)|, |??"(x)|} < b 2 . 

2. The initial state satisfies (v a ,p a ) £ Hq cr (f2;M JV ) x fl where 

K := ju £ R 1 (n) \ ipi < v < 1)2 a.e. in flj , 

with — 1 — TpZ =: ip! < 0 < ip 2 ■= 1 — 1 Ta- 

3. The density p depends on the order parameter p via 

p(p) = max | P 1 + P 2 + P 2 ^ Pl (v? + ^),o| > 0. 

We note that by Remark 2.1 the pure phases are attained at x when p(x) = if 1 
or p(x) = ip 2 , and the max-operator in Assumption 2.2.3 ensures that the density 
remains always non-negative. The latter is necessary to derive appropriate energy 
estimates. 

With these assumptions we now state the semi-discrete Cahn-Hilliard-Navier- 
Stokes system. For the sake of generality, we additionally introduce a distributed 
force on the right-hand side of the Navier-Stokes equation, which will later serve the 
purpose of a distributed control. Below and throughout the paper, r > 0 denotes 
the time step-size and M £ N the total number of time instances in the semi-discrete 
setting. 

Definition 2.1 (Semi-discrete CHNS-system). Let T 0 : H (f2) —> R be a convex 
functional with subdifferential 9’IV Fixing (</?_i,uo) = (p a ,v a ) we say that a triple 

(w,«) = (fe)f.o 1 ,(ft£ 1 ,(€( 1 ) 
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in x H 2 dri (n) M x i?Q R Ar ) M_1 solves the semi-discrete CHNS system 

with respect to a given control u = (ui)^.] 1 £ T 2 (f2;R Ar ) M ~ 1 , denoted as (ip,p,v) £ 
Sq,(u), if it holds for all (f> £ H 1 (Q) and if £ Hq a (It; l w ) that 


( 2 . 1 ) 

( 2 . 2 ) 


(2.3) 


^ —^ + ( Wi+1 v^ i5 0) - (div(m(</?j)V/ij+i), 0) = 0, 

(—Ayjj+i, </>) + (d^i+i)A) - {m+iA) - {wiA) = o, 

/ P(<Pi)Vi+l - p(<Pi-l)Vi \ , , , . , 

< -> V ) + (div(u i+ i <g> p(tpi-i)Vi),ip) H -i H 1 

\ T / Hr 1 .Hi °’ <T °’ a 


- ^div(u l+ i 


P2 - Pi 




H - 1 ,Hi, 

0. cr ’ 0 ,cr 


(2ri(tpi)e(vi +1 ),e(il>)) 


(pi+l^Ti, ' 4 ’)H - 1 H 1 ~ ( u i+li 4 ’) H " 1 .H 1 

U ,<J ' U, cr U,ct' U,( 


TTie first two equations are supposed to hold for every 0<i + l<M — 1 and the last 
equation holds for every l<i + l<M — 1. 

Remark 2.3. /n general, the subdifferential of a convex function d'o can be 
a set-valued mapping, see, e.g., [15]. In this case, by equation (2.2) there exists 
/3 £ c?\I/o (v^i+i) such that 

(-A (p i+1 ,<j>) + - {pm,(f>) - (mPiA) = 0 , £ h 1 ^). 


We note that in the above system the boundary conditions specified in (1.1) are 
included in the respective function spaces. 

It is interesting to note that our semi-discretization of (1.1) in time involves three 
time instances (i — l,i,i + 1). Equations (2.1) and (2.2), however, do not involve the 
velocity at the ’’old” time instance i — 1. As a consequence, (po, pf) are characterized 
by the (decoupled) Cahn-Hilliard system only. At the final time instance, however, the 
coupling of the Cahn-Hilliard and the Navier-Stokes system is maintained; otherwise, 
we have little hope to derive some energy estimates for the system. 

Finally, we present the optimal control problem for the semi-discrete CHNS sys¬ 
tem. For its formulation, let U a d C L 2 (H;M Ar ) M_1 and J : X —> R be a Frechet 
differentiable function, with 

X :=h\ n) M x H\n) M x Hq a (fl; M jv ) m_1 x L 2 (fl- R^" 1 . 

Further requirements on U a d and J are made explicit in connection with the existence 
result, Theorem 4.1, below. 

Definition 2.2. The optimal control problem is given by 

^, min J(tp, p,v,u) over (p, p,v,u) £ A 

s.t. u £ U ad , (AMA) e Sm(u). 

In many applications, J is given by a tracking-type functional and U a d by uni¬ 
lateral or bilateral box constraints. 
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3. Existence of feasible points. In this section, we prove the existence of 
feasible points for the optimization problem (P^). As stated earlier, for deriving 
stationarity conditions we will later on approximate the double-obstacle potential by 
a sequence of smooth potentials of double-well type. Therefore, we consider here the 
following two types of free energy densities. 

Assumption 3.1. The functional To : H 1 (f2) — > R is convex, proper and lower- 
semicontinuous. It has one of the two subsequent properties: 

1. Either it is given by To(<^) := ipo(tp(x))dx where ip o : R —>■ R := 111 {+oo} 
represents the double-obstacle potential, 


( +00 

if z <ipi, 

iPo(z) := < 0 

<N 

i5- 

VI 

0 

VI 

[ +00 

if z > ip 2 . 


2. Or it originates from a double-well type potential and satisfies: 

(a) To is Frechet differentiable with (Tq(</?)} = <9To(</?) C L 2 ( SI) for every 

V £ H 1 { fi); 

(b) There exists B u £ R such that To(v?) < B u for every tp £ K. 
Additionally, we assume that the functional T(tp) := To (p) — f n ^<p(x) 2 dx, n > 0, 

is bounded from below by a constant P; £ R. 

We start by studying the semi-discrete CHNS system for a single time step. For 
this purpose, assume that the pair (ip,v) £ H (SI) x Pq j(T (S2; R N ) is given. We then 
show the existence of a point (<p, p, v) which solves a slightly modified system. The¬ 
orem 3.5 collects the results for all time steps via an induction argument. Finally, 
Theorem 3.7 shows that the modified system equals the original CHNS system under 
suitable assumptions. 

The starting point for our considerations is an energy estimate for the generalized 
system. This estimate will be useful to establish the boundedness of the feasible set. 
We note that in what follows, C, C\ and C 2 denote generic constants which may take 
different values at different occasions. 

Lemma 3.1 (Energy estimate for a single time step). Let dp £ H (SI), v £ 

0„ £ 0„,0 V £ H~\n), zz £ 1 £ 

P 2 (S1), /o,/_i > 0 be given such that 


(3.1) 


fo - f -1 


+ divzz = 0 a.e. on SI. 


In case of the double-obstacle potential suppose additionally that dp £ K. 
Then, if (ip,p,v) £ H (SI) x H (S2) x ifg_ cr (S2;R JV ) solves the system 


(3.2) 

(3.3) 


^ t ^ ,<t>J + {vV<P,<P) ~ (di v(m(<p)Vp),(p) = (0 M ,<£), V0 £ 

- (/b <P) - {k<p, <p) + (—Ay>, <p) + (< 9 T 0 {p),(p) = ( <d v , <p ), £ H l (Vl), 

— -+ (div{v®v),ip) H -i „i +(2r)(<p)e(v),e(ip)) 

T / h~1,h 

— (pV(p,ip) H -i H i = (Q v ,ip) H -i H i , V'0 £ H 0 (J (Sl; R ), 

U.cr ’ U ,(7 U ,g ’ U ,T 


(3.4) 
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the following energy estimate holds true: 


fo M' 


■dx ■ 


I 


■dx + ^(<p) + / f—i 


\v — V 


2 dX + Jn 


IVip-Vtf 


dx 


+ t / 2r](ip) \e(v)\ 2 dx + t / m(tp) \V p \ 2 dx + / k 


(v - vf 


(3.5) 




iv#r 


dx + ^(p) + g(p,p,v), 


where g is defined as 


(3.6) 9{v,li,v) := (0^,/z) + ( 0 y , ^ ^ ^ 'j + (0«,w)jj-i iff i ( 


Proof. First, we observe that 

(div(u ® i/), v) = ((divi/)u + (i/ • V)u, v) 

f . v f v 

= / ((divi/)— + (y ■ S/)v)vdx + / (div^)—uda: 

J n 2 J n 2 

= [ div (#) + (divO^dx = [ (drvuKdx. 
Jn \ 2. J 2 J n 2 


(3.7) 

Next, one verifies 

(fov- f-iv,v) = 

(3.8) 


/o M 


+ 


/-i M' 


■dx — 


(fo ~ f- i) M 


■da: 


■da; 




■da;. 


Testing (3.2),(3.3) and (3.4) with /a, ^ and v, respectively, summing up and inte¬ 

grating by parts, we obtain 


0 = 


fo M 2 - /-i l^l 2 

2 r 


da: + / 

an 2r 


■dx 


(fo - /~i) kr 

2 r 


da; 


(3.9) 


+ f (divzd)-^j-da: + f 2ri(ip) \e(v)\ 2 dx + f m(ip) |V g\ 2 dx 
Jo, 2 Jq Jq 

1 f (p — (p 

+ - (d'k 0 (ip), (f - (jfi-jj-i jji - k / -da; 

T ’ J n T 

+ ~ Vip(Vip -Vtp)dx - g(ip,n,v), 

T Jn 


where we also use the previous equations (3.7) and (3.8). From the definition of the 
subdifferential we infer 

(3.10) (d^ 0 (p), ip - <p) > T(^) - 'F(^) + ~ f^dx. 

Inserting (3.1),(3.10) into (3.9) and using 2 a (a — b) = a 2 — b 2 + (a — b ) 2 once for 
(a, b ) = (Vp, V<)5) and then for (a, b) = (ip, ip) we obtain the assertion. 0 
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Remark 3.2. Note that the system (3.2)-(3.f) corresponds to the system (2.1)- 
(2.3) for one time step only when choosing 

v = v it <p = <Pi, fo = p(<Pi), f-i = p{<Pi-i), 
v = p(<pi-i)vi - P2 Pl pi, 

QLp - - O' 

Now we prove the existence of solutions to the system (3.2)-(3.4). The proof mainly 
relies on the application of Schaefer’s fixed point theorem, also called the Leray- 
Schauder principle, and combines arguments from [1, Lemma 4.3] and monotone op¬ 
erator theory. 

Theorem 3.2 (Existence of solutions to the CHNS system for a single time 
step). Let the assumptions of Lemma 3.1 be satisfied. Then the system (3.2)-(3.f) 
has a solution (ip 7 p,v) £ .ff^fl) x R 1 (f2) x Hq CT (fi; S. N ). 

Proof. We start by defining 

(3.If) X := H\n) x h\lI) x R N ), 

(3.12) Y :=H~\n) xH~\Ll) x R N )*, 

and the operators Q\ : H 1 (Ll) —► H 1 (Ll), Q 2 : H 1 (Ll) H 1 (fl), Q 3 : Hq a (Ll; R w ) —>■ 
Hq a (Ll; 1^)*, 0 : X =4 Y and LF : X —>• Y (here and below ’=£’ indicates a set-valued 
mapping) via 

Gi(p) ■■= -di v(m(ip)Vp) - 0^, G2{<p) : = + d'&oiv) - 0 V , 

Gz(v) ■= -div(2ry(^)e(i;)) - 0„, 

G(y,p,v) := (Gi(p) 1 G2(t),G3(v)) T , F(p,p,v) := (J*i, F 2 , T 3 ) t , 

with 

Ti(p,tJ,,v) := --- -vS7<p, T 2 (ip,p,v) := p + K<p, 

r 

J 3 O p, £t, v) := ——- Lj— — div(u <g» v) + pVp. 

T 

Using this notation, the system (3.2)-(3.4) can be stated as 

(3.13) 0 £ G{<p, v) — LF(<p, p, v) C Y. 

By standard arguments, the mappings Gi and Q 3 are invertible and the respective 
inverse mapping is continuous. Since the Laplace operator is invertible from H (f2) 
to H (LI) and the subdifferential ch&o is maximal monotone (cf. [38, Theorem A]), 
Q 2 is invertible, as well. Concerning the continuity of Gf 1 , let £ 1,^2 £ H (Ll) and 
<Pi, tp 2 £ H (Ll) satisfy ipj = G 2 l (fj) for j = 1, 2. Using Poincare’s inequality and the 
monotonicity of <9To, we immediately obtain 

11^2 - <Pi\\ 2 h i < C((-A(ip 2 - <Pi),<P2 - <pi) + (d'&oiw) ^ d^ 0 (pi),ip 2 - <^i)) 

— C (f 2 — (i,ip 2 — ipi) < C ||^2 “ CiIIjy —1 11^2 — <Al||#i > 

showing the continuity of Gf 1 ■ 
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Due to the compact embedding of the space Y := L,i(fl) x L%(il) x (Cl-,M. N ), 
into Y, the inverse of Q is a compact operator from Y to X. Further, T : X —> Y is 
continuous. Hence, the operator T o Q~ 1 : Y —> Y is compact. 

In what follows, we show the existence of a solution 6* to the fixed point equation 

(3.14) 6* - Tog- 1 {8*) = 0 £ Y. 

Then it immediately follows that Cfyfyd*) solves the system (3.2)-(3.4). In order to 
apply Schaefer’s theorem with respect to the operator J- o Q we verify the condition 
that the set D := (Jo<A<i e = A.F ° £ _1 (^)} is bounded. For this purpose, 
assume that 5 £Y and A € [0,1] satisfy 

(3.15) 6 = \Tog~ 1 (6), 

and define (ip,fi,v) := Q~ 1 (S) £ X. Thus, (3.15) can be rewritten as 

(3.16) g(ip,fi,v) - XF(tp,n,v) = 0 


which is equivalent to the following system of equations 


^ ^ ,<t>j + (AvVcfy 4>) = (div(m(<^)V/i), <j>) + (0^, </>), V</> £ iT^fi), 

<A^t, 0) + (\Ktf,<j>) = (-A ip,(t>) + (d’Fo (v)),4>) - (©*,,<£), \fy> £ H 1 ^), 


a /——+ A (divfy <g> v), 
\ T / h~Ih^ 


= \(fj,V<p,i/))H-i h' +( q v,^)h- 1 h 1 , £ H 0 (Cl; R ). 

0,cr ’ 0, cr 0,<7 ’ 0, cr 

Analogously to the proof of Lemma 3.1, we test this system by /i, and u, respec¬ 
tively, sum up the resulting equations and integrate by parts to derive 


0 = A 

(3.17) 

which leads to 


h I”' 2 - + A [ 


2 T 


/ n 


2r 


2ry(^) |e(u)| da: 


+ [ m{(p) \X fi\ 2 dx -\— [ d^>o(ip)(if — tp)dx — \n / f———dx 

Jn T Jn Jn T 


+ - [ Xf(Xip - V<f)dx - g(ip, f, v), 

t Jn 


(3.18) 


2rj(if) \e(v)\ 2 dx + [ m(<p) \ Xg\ 2 dx + -'F(fy) + - I \S7f\ 2 dx-g(if,g,v) 

! Jn t t Jn 

If 0 1 

dx H— / |V<p dx-\ — ^(f)- 

T Jn T 


f -1 |i?r 

2 T 


Note that for obtaining (3.17) we also make use of (3.1). The right-hand side of (3.18) 
can be bounded by a constant C := C(N, H, r, f)> 0 which is independent of 

A. Since ’F is bounded from below, this leads to 

(3.19) j 2r)(if) \e(v)\ 2 dx + f m(if) | V fi\ 2 dx H— j |V(/?| 2 dx < C + g(if, /a, v). 

J Q, J Q J Q, 
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Due to Korn’s inequality, Poincare’s inequality and from the boundedness of r)(-) and 
m(-), we infer 

Mm + \\pfm + Mlm < C + g(y,p,v) 

(3.20) < Ci + + UmIIhi + 

where C 2 > 0 depends only on 0 M , 0 V and 0„. The last inequality yields the 
boundedness of (<p, p,v) in X. Next, we derive bounds for T. In fact, we have 

CM + IHI + II£M, 

||J r 2 (v5,/i,r;)|| i 3/2 < C(||/i|| + II^H), 

\\X 3 {ip,ii,v)\\l3/2 < C(||u|| ffl + ||u|| H1 IMIffi + INI M\m + INIhO- 

Since (p , v and v are fixed, D is bounded in Y. Hence Schaefer’s theorem is applicable 
implying that equation (3.14) admits a fixed point 6* £ Y. Then G~ 1 (S*) solves the 
system (3.2)-(3.4). □ 

In our setting, the right-hand sides of the system (3.2)-(3.4) are square integrable 
functions. This enables the derivation of higher regularity properties for the solutions 
obtained in Theorem 3.2. 

Lemma 3.3 (Regularity of solutions). Let the assumptions of Lemma 3.1 be 
satisfied, and suppose additionally that 0 M ,0 V £ L 2 {Cl), as well as /o,/-i £ L 3 (Cl) 
and (p £ H 2 (Cl). 

Then it holds that ip, p £ Hg n (Cl) and v £ H 2 (Cl','E. N ), provided that (ip,/j,,v) £ 

H (f2) x H (Cl) x Hq a (Cl; M. N ) satisfies the system (3.2)-(3-4). Moreover, there exists 
a constant C = C(N, Cl, b±, & 2 , t, k) > 0 such that 

IMIif2 + IImII h 2 + IMI,?/ 2 

( 3 . 21 ) < C(M + ||/i|| + |H| + || 0 J + ||©pt|| + Mm Mm + K(v>)ll)- 

In case of the double-obstacle potential, it also holds that p £ K and the term ||4/g(<^)|| 
in the above inequality is dropped. 

Proof. Equation (3.3) is equivalent to 

(3.22) Atp + gi £ d'Yo(p) in H 1 (D) 

with gi := p + ntp + By Sobolev‘s embedding theorem g\ is in L 2 (Cl). In case of 
the double-well type potential, Assumption 3.1.2 (a) then implies g 2 := — gi + 'h'^p) £ 

L 2 (Cl) and Aip = p 2 . Applying [34, Theorem 2.3.6] and [34, Remark 2.3.7], we deduce 

_2 

that ip £ H (Cl) is the unique solution of the Neumann problem 

Aip = g 2 in Cl, d n <P\dn = 0 on dCl. 

Furthermore, [34, Theorem 2.3.1] yields the existence of a constant C := C(N,Cl) 
such that 

(3.23) |M|* a < C(M + || ff2 ||) < C(M + ||m|| + « |H| + ||©„|| + \\%(P)\\). 

In case of the double-obstacle potential, (3.22) is equivalent to the variational inequal¬ 
ity problem: 


(3.24) 


Find ip £ K : (—Aip — g\, <f> — ip) > 0, V</> £ K. 
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Then the assertion follows from the subsequent lemma. 

Lemma 3.4. If p £ K solves the variational inequality problem (3.24) with g± £ 
L 2 (ft), then p £ Hg (ft) and there exists a constant C = C(N,Q) > 0 such that 
\\<f\\ H 2 < CUsili . 

For the sake of completeness we provide a proof in the appendix. It closely follows 
the lines of argumentation of [31, Chapter IV]. 

Regarding p, we argue similarly. Indeed, first note that by Sobolev's embedding 
theorem and Holder’s inequality <73 := + vVp — — p is an element of L 2 (ft). 

Furthermore, the coefficient function m(p) is contained in 7F 2 (ft) and IV 1 , 6 (ft), re¬ 
spectively (cf. [31, II, Lemma A.3]). Equation (3.2) is equivalent to 

(3.25) m{p)Ap + V(m(y5))V/i — p = <73 in H (12). 

Again by [34, Theorem 2.3.5] and [34, Theorem 2.3.1] p £ H^ n (Cl) and it holds that 

(3.26) \\p\\ H 2 < C(\\p\\ + || 5 3||) < C(\\p\\ + |M| + ||p[| + |M | h1 \\p\\ H2 + ||© M ||), 

where C > 0 depends only on TV, ft, 61 , 62 , r. 

Finally, we show the desired regularity of v. Since div(e(w)) = i(Aw + V(divw)), 
it holds that 


div(2ri(p)e(v)) = 2p' (p)e(v)Vp + rj(p)Av, 
and therefore by equation (3.4) that 

(3.27) Aw = rj{p)~ l div(w ® v) — 2rj'(p)e(v)V p — (/qw — f-iv) — pVp — O t 


Moreover, div(w 0 v) = (Dv)v + wdivjA By the assumptions all summands in the 
second line of (3.27) belong to L 2 (f2;R JV ) and v £ II 1 (fl; R^). Hence, we have 


(3.28) Av = r](p) 1 [{Dv)v + vdrvv — 2rj'(p)e(v)Vp + /], 

with Il/H < C(~) for a constant C(z) > 0 depending only on 


2 = (A, ft, 77 , t, \\p\\m, ||<^||jf 2 , \\p\\ H ^ ||0«||). 

In order to show v £ iF 2 (ft;R JV ), we apply a bootstrap argument and well-known 
regularity results for the stationary Stokes’ equation, cf. [41], 

1. Since v £ Hl a (VL\R N ), we have that ( Dv)v , e{v)\7p and vdivv belong to 
L 3 / 2 (ft; K^). Therefore, [41, Prop. 2.3, p. 35] and (3.28) show that v £ IF 2 ’ 3 / 2 (ft; R^) 
and that ||w|| W 2 , 3/2 < C(z ) for a constant C depending only on z. 

2. Next, v £ kF 2 , 3 / 2 (ft;R^) and the continuous embedding of VF 1 , 3 / 2 (ft) into 

L 3 (ft) (which we denote by kF 1 , 3 / 2 (ft) L 3 (ft)) imply that (Dv)v and e(w)V<p 

belong to L 2 (ft;R Ar ). Moreover, IF 2 ’ 3 / 2 (ft) L p (fl) for every p < 00. Hence 

vdiw £ L 2 _e (ft;R Ar ) for every e > 0. Applying [41, Prop. 2.3, p. 35] again yields 
v £ IF 2 ’ 2 -£ (ft; R w ) for all e > 0 and ||w||t^ 2 , 2 - e < C(e, z). 

3. Finally, having v £ VF 2 , 2 _e (ft;R w ) and since IF 2 ’ 2 -e (ft) °->- L°°(ft) for e 
sufficiently small, it follows that also vdivv belongs to L 2 (ft;R Ar ). Thus, we arrive 
at v £ iJ 2 (ft;K Ar ) and ||w ||#2 < C(z). 

This completes the proof. □ 
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Our aim is to prove the existence of a solution to the semi-discrete CHNS system 
with the help of Theorem 3.2. This result, however, is not directly applicable with the 
setting of Remark 3.2, as /o, /_i and v need not satisfy equation (3.1). This is due to 
the nonsmoothness of the density function and the fact that ip may attain values in 
R. (rather than [ipi,^]) for a double-well-type potential. We overcome this difficulty 
by applying Theorem 3.2 with the setting 

^ v := Vi , if := fi, f 0 := p(fi), /_i := p(fi-i), 

v := V{vi , ifi, fi- 1 , pi), <d v := u i+1 , := 0^ := 0, 

where v : H)j a (Q; R w ) x H 2 (Ll) 3 -*• H 1 {Q-,R N ) is given by 


(3.30) v{v,f,(p,p):= 

Here G : L 2 (Ll) 


p{f)v — P2 2 Pl m(<p)V p if p(p), p(f) > 0 a.e. in fi, 
G (£M^pM) else. 

i i —> Ci is a solution operator to 


(3.31) —div<C = 5 a.e. on Ll. 

A specific realization of G is obtained by first solving — = 6 in L 2 (Ll) with £ = 0 

on dLl, yielding £ £ H 2 (LI) (~l Hq(LI), and then setting ( '■= £ H 1 (Ll,R N ). We 

next prove that there always exists a solution to the system (2.1),(2.2),(3.32) where 
the semi-discrete Navier-Stokes equation (2.3) is replaced by 


(3.32) 


p(lfi)v i+ i - p(fi-i)Vi 




' T 7 
+ (div(v i+ i fi, fi-!, p^),ip) H -i H i - 

0,(7 ’ 0 , <7 

= {ui+i, $)h-^,hI „ 


- (2ri(fi)e(v i+1 ),e{i/j)) 

(Pi+l^fii ‘4’)HZ 1 ,H 1 

0,(7 ’ 0 , CT 

, vv-£i?oV( fi ; Rjv )- 


We point out that (3.32) coincides with (2.3) if 


(3.33) 


min{p(^), p(pi_i)} > 0 a.e. on LI. 


For the double-obstacle potential this always holds true, since ft is contained in the 
interval Hence in this case p(fi) > p(if i) = pi > 0. 

In Theorem 3.7 we show that for the double-well type potentials under consider¬ 
ation f remains in a close neighborhood of [ifi , 1 P 2 ] and therefore condition (3.33) is 
satisfied as well. 

Proposition 3.5 (Existence of solution to a modified state system). Let V : 
Hq CT (fl; R^) x H (LI) 3 —> i7 1 (H;R Ar ) be defined by (3.30). Then for every u £ 

L 2 (Lt'M. N ) M ~ l there exists apoint (ip,p,v) £ Hg n (Ll) M xHg n (Ll) M xHq cr (H;R Ar ) M_1 
which solves the semi-discrete system (2.1),(2.2),(3.32). Moreover, every solution 
(f,p,v) satisfies (f,p,v) £ H 2 g n (LL) M x i7g n (H) M x H^^Ll; R^^ -1 . 

Proof. Standard arguments guarantee the existence of (p 0 ,p 0 ) £ H (H) x H (LI) 

_2 

such that (2.1)-(2.2) is satisfied for i = —1. Lemma 3.3 yields (ipo,po) £ H dn (Ll) x 

H 2 aM- 

If condition (3.33) holds true, then Assumption 2.2.3 and (2.1) imply 
di vV(vi,fi,fi-i, pi) = P2 Pl (Wfi-iVi - dYv(m(ifi-i)Vpi)) 

= P2 ~ Pl . fi ~fi-1 . = p(fi) - P(fi-l) 

2 T T 
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Consequently, if (ipt, Pi,vf) satisfies (2.1), then assumption (3.1) is always satisfied by 
the definition of V in the sense that 

(3.34) Pifi) - P&i i) -)_ divX7(^, tpi, Hi) =0 a.e. on fh 

T 

Therefore, we can apply Theorem 3.2 with the setting (3.29) for i = 0 to guarantee 
the existence of (ipi, /ii, iq) £ H (f2) x H (Q) x Hq ^(fl; R w ) such that the system 
(2.1),(2.2),(3.32) is satisfied. 

Now Lemma 3.3 yields x H^ n (Q) x In the 

case of the double-obstacle potential it additionally follows that ipi £ K. 

Repeated applications of Theorem 3.2 and Lemma 3.3 for each time step i = 
1,.., M — 2 prove the assertion!] 

As discussed above, this theorem directly guarantees the existence of a solution 
to the semi-discrete CHNS system for the double-obstacle potential. Next we address 
the boundedness of the solutions which is needed later on to ensure the existence of 
optimal points for (P$). For this purpose, we apply the energy estimate of Lemma 
3.1 at each time step. 

Lemma 3.6 (Boundedness of the state). There exists a positive constant C = 
C(N, S2, bi, 62 , t, k, v a , u) > 0 such that for every solution (<p,/j,,v) £ (Q) M x 
H dn (n) M x Pq ^.(11; R jv ) m_1 of Theorem 3.5 it holds that 

(3.35) IMI(ff 2 )M + ||M||(J/ 2 )M + ||<p||( ff 2 )M + l < C. 

Furthermore, the operator L 2 (tt, K N ) M_1 9 u i —> C(N, fl, b\, 62 , t, k, v a , 'pa, u) £ K 
is bounded. 

Proof. We define the functional E : Hq a (£l) x H 1 (f2) x H 1 ( f2) —> K. as follows: 

(3.36) E(v, <t>) := ( dx + f dx + ^/{p). 

J n z Jn 1 

Let j £ {1, ,.,M — 2} be arbitrarily fixed. Then by repeatedly applying Lemma 3.1 
with the setting (3.29) for i = j,j — 1,.., 0, we conclude that 

E (vj+i,Vj+i,Vj)+ r 2p(ipj) \e{v j+ i)\ 2 dx + t m(tpj) |V p j+1 \ 2 dx 
J n Jn 

< E{vj,ipj,^j- 1 ) + {u j+1 ,v j+1 ) 

< E(Vj-i, tpj-!, ipj- 2 ) + (Uj,Vj) + (u j+ i,v j+1 ) 

3 +1 

< E(v 0 , <po, tp- 1 ) + ^2 ( u £ Vi ) • 

i= 1 

By Assumptions 2.2 and 3.1 this yields 

/ ^ (f3+1 ^ dx + n) + 2rbi [ \e{vj +1 )\ 2 dx + rbi [ |V/Xj+i| 2 dx 
Jn 1 Jn Jn 

M -1 

(3.37) < E(v 0 ,<p 0 ,<p-i) + ll“*ll INI < Ci + C 2 ||u|| (i 2 ) M-i , 

i=1 






M. Hintermiiller, T. Keil, D. Wegner 


15 


where C\ depends only on the initial data (N,Q, Bi, B u ,v a , ip a ). Due to Korn’s in¬ 
equality and Poincare’s inequality, this ensures 

\\v j+ i\\ H 1 + llAb'+l II// 1 + ll^’j+lllfl-i < Ci + C 2 \\u\\ ( L 2 )M-i |M|(i/l)M ■ 

Since j £ {1 ,M — 1} is arbitrarily chosen, we infer 


(3.38) IMI(jji)M + IImII(hi)" + IMI(/r 1 )*f+ 1 < Ci + C 2 ^||w||( L 2)m-i IMI(//i)m) ■ 

Hence (ip, /x, v) is bounded in H 1 {VL) m xH 1 (Q) M xHq^^LI; R w ) m_1 . Then the bound¬ 
edness in the respective U 2 -spaces follows directly by applying Lemma 3.3 for each 
time step.D 

Next we address the case of the double-well potential and show that for appro¬ 
priate double-well type potentials the order parameter of a solution to the system 
(2.1),(2.2),(3.32) is always greater than ipi — e for some small e > 0. 

Theorem 3.7. Let u £ L 2 (fl; R w ) u_1 be given and a sequence of 

l J k(E N 

functions which satisfies the following two conditions: 

(k) 

1. For every k € N fulfills Assumption 3.1. 

2. If N is a sequence in H 1 (fi) such that there exists C > 0 with 

< C for k £ N, then 


nm(-^^ +^i,0) 


L 1 


0 , as k —> oo. 


Furthermore, let w^ fe ^)} fcgN be a sequence of solutions to the systems 

(2.1),(2.2),(3.32) with ^ ■ Then 


max(-^ (fe) +^i,0) 


0 , as k —> oo. 


Proof. Employing Lemma 3.6, in particular inequality (3.37) from its proof, then 
we see that for every i £ {—1, ..,M — 1} and k £ N it holds that ^^(ip^) < C\. 
Hence, we conclude 




(fe) 




<Ci + 




(k) 


L2 


< c 2 . 


By assumption, this yields 


(3.39) 


max(-^ fc) + ipi, 0) 


L 1 


—>• 0 for k —>• oo. 


Next, we use the technique of [24, Proposition 2.4] and [24, Remark 2.5] to derive 
that max(— tp^ + ip i,0) —> 0 for k —> oo. We stay brief here and refer to [24] 

for details on the technique. By Lemma 3.6 the sequence N is bounded in 


—2 


H (fi), and due to Sobolev’s embedding theorem it is also bounded in W 1,6 (H) and 
C 0,/3 (H), (3 < \, respectively. Thus, there exists a constant Cp such that for every 
k £ N we have ||v?^ || COi/3 < C/3- 
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For fix k £ N assume that 


max(-pf ) +^i,0) 


> 0 and define the set G := 


jw £ : (p^ k \u) <if >l < oj. Then let u) max £ G satisfy 


-<Pi k \uma X ) + V'l — -<Pi K> + 4>1 


(fc) 


L°°(G) 




l~( n) 


Due to the Holder continuity of tp\ k \ for every x £ H which satisfies |x — w maa: | R Ar < 


2C« 


it holds that 


-ip\ k) (x) +tpl> -tpf’ (Umax) + V’l “ 

S k )( 


( fc ) / 




(fc) 


C(0./3)(O) 




Vi (Umax ) + 'fil > g 


As fl satisfies the cone condition, there exists a finite cone K r (uj ma x) '■= K(u max ) (~l 
B (tUmax, r) of radius r and with vertex u max such that K r (u ma x ) C LI. Hence the cone 


Kr(uj max) with R := min ^r, ^ Vi ) I is contained in G. Consequently, 

we find 


max(-<£- fc) + ^i,0) 


> 


k (^) JK R (ui max ) 


> 


> 


+ tpidx 

(-^'(Wmrn) + V’l) 


dx 


K R (jJ^max ) 


max(-(jf^ + V’l, 0) 


L°°(n) 


In combination with (3.39) this proves the assertion. □ 
We define ip~ £ R as 


(3.40) 


tp := inf {tp £ R : p(ip) > 0} < V’l- 


Let u £ L 2 (Ll; R JV ) M 1 be given and let { 1 be a sequence of double-well type 

potentials which approximates the double-obstacle potential in a certain sense, i.e., 
it satisfies condition 2 of Theorem 3.7. Then Theorem 3.7 ensures that there exists 
k* £ N such that for every k > k* the solutions (< p^ k \p^ k \v^) to the corresponding 
systems (2.1),(2.2),(3.32) with T 0 = ’Fq^ satisfy 

(3.41) ip^ >¥>", Vi = -1,..,M-1. 

Hence p(Vi^) > 0 for every i = —1, ..,M — 1 and k > k* . Thus, (3.32) coincides with 

(2.3), which leads to the subsequent theorem. 

Theorem 3.8 (Existence of feasible points). Let u £ L 2 (H;K JV ) M_1 . Let 4/o 

be the double-obstacle potential defined in Assumption 3.1.1 and let < \ be a 

l J fceN 

sequence which satisfies the conditions of Theorem 3. 7. 
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Then there exists k* £ N such that the system (2.1)-(2.3) admits a solution 
fp,p,v) for every 4'o £ {^0} U j . For this solution (ip,p,v) the result 

of Lemma 3.6 remains true. 

In other words, the semi-discrete CHNS system (2.1)-(2.3) has a solution if the 
double-well type potential under consideration is close enough to the double-obstacle 
potential. In the following sections we always assume that this is the case. In Defi¬ 
nition 7.1 below, we propose a specific regularization which satisfies the conditions of 
Theorem 3.7. 

4. Existence of globally optimal points. The previous section guarantees the 
existence of feasible points for the optimal control problem (P^). Next we investigate 
the existence of a solution to (Pip). For this purpose, we need to impose additional 
assumptions on the objective functional and the constraint set U a d- 

Theorem 4.1 (Existence of global solutions). Suppose that J : H 2 g n (Vt) M x 

H 2 g n (Ll) M x Hq a (Q- 1 K Ar ) M_1 x L 2 (fl;M Ar ) M_1 —>• R is convex and weakly lower-semi- 
continuous and U a d is non-empty, closed and convex. Assume that either U a d is 
bounded or J is partially coercive, i.e. for every sequence p^ k \ v^ k \ u^)} fceN 

with limfe^oo ||«(*) || = 00 it holds that lmp.-^ p^ k \ v^ k \ u^) = 00 . Then the 

optimization problem (Pq,) admits a global solution. 

Proof. By Theorem 3.5 the feasible set of the problem (Pp) is non-empty and 
contained in H 2 gjn) M x Pa n (D) M x Hq R N ) M ~ 1 x U ad . 

Let {(tp( k \ p( k \ y( k \ u ^)} ken be an infimizing sequence of J in x 

H 2 dn (n) M x R Ar ) M_1 x U ad with £ ^(#1) such that 

(4.1) lim J(<p< k \pW,vW,u w ) = inf J(ip,p,v,u). 

fc-> 00 ueu ad ,(<p,ij.,v)&s^(u) ' 


Note that the infirnum on the right-hand side may be — 00 . The sequence {u( fc )} fceN 
is bounded in the reflexive Banach space L 2 (f2;M JV ) Af . This follows either directly 
from the boundedness of the set U ad or from the partial coercivity of J. Then 
by Lemma 3.6 the sequence {^ k \ p^ k \ v^) is bounded in H 2 g n (Q) M x H 2 g n (Tf) M x 
Hq a (Q,] R Ar ) M_1 . Setting {w (/c) } fcGN := {{^ k \ p^ k \v ^, u (fc) )} fcgN , there" exists a 
weakly convergent subsequence with limit point w* := (ip*, p*, v*, u*) £ 

H 2 g n (f V) M x H 2 g^ (f T) M x Hq CT (f l; R JV ) M_1 . Using the weak lower-semicontinuity of J , 
this implies 

—00 < J{w*) < liminf J(w^ kl ' > ) = inf J(ip,p,v,u) 

l —>00 u£U a di( l P,L L ) v )£S'i'(u) 


where the last equality holds due to (4.1). Since U ad is weakly closed, u* belongs to 

Uad- 

It remains to show that (ip*, p*,v*) £ Sy (it*). For this purpose, we write l instead 
of ki, and we start by considering the limit of div(u -^ 1 ® P2 ~ P1 m(p^ 1 )'V p^), 
for arbitrary i £ {0,.., M — 2} and if £ P 1 (D; R^). Using the triangle and Holder’s 
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inequality we derive 

m(^l i)V/ 4 0 • Vi/j - m(<p*_ 1 )'Vn* ■ Vtf 


- v/<) • 


L 4 /3 

,(0 


i(d-i) 


Since V/i -^ converges weakly to V/x* in H 1 ^) and iJ^fl) is compactly embedded into 


V/4°- V/x* 


L 4 


L 4 /3 

l|V^|| i2 


( m Oi-i) - TO (<P*-i))V/x* ■ Vx/> 


m (V 5 i-l) — TO (^i-l) 


L 4 /3 

ll v M*IL 4 l|VV’|| i 2 • 


L 4 (H), 


V/xf > - V/x* 


L 4 


tends to zero for l —>• oo. Due to the compact embedding of 


H 2 (Tl) into W 1 ’ A (Cl), we have —$■ P*-\ strongly in kF 1 , 4 (f2). Due to Assumption 
2.2.1, m is Lipschitz continuous. Since W 1,a (VL) can be embedded into L°°(f2), we 


infer 




0 . 


Consequently, the sequence p2 2 Pl m(<p[ l \)'S/■'V i/j converges strongly in L 4 / 3 (fi) 
to P2 ^ Pl m(p*_ 1 )\' p* • Vx/>. By Sobolev’s embedding theorem and the weak con¬ 
tinuity of the embedding operator, nSj converges weakly in L 4 (f2) to v* +1 . Hence 
^div(u i ^ 1 0 p2 ~ pi m(y?^ 1 )V/x^), -0^ converges to (div(n * +1 0 P2 ~ Pl w(yj*_ 1 )V/x*), 
as l —>• oo. 

One proceeds analogously for the remaining terms in the system (2.1)-(2.3) which 
do not depend on the subdifferential of To- 

In this way we also show that A + pl+i + k ipf'* converges strongly in H (f2) 
to A (f* + i + p* +1 + wp* for every i = —1,.., M — 2. Furthermore, Pi+i —> p* + i in 

H l (Vl), and for every l G N it holds that A ip^ + p^i + Kipf^ € S^o^i+r)- Due to 
the maximal monotonicity of cITo, this implies 


(4.2) Ap* +1 + p* +1 + Kip* G 3*o(v>* +1 ) 

for every * = — 1,.., M — 2. In summary, we have shown (ip*, p* ,v*) G S\s/(u*). Hence 
the w* is contained in the feasible set of the problem (Pj) and therefore solves the 
problem!] 

5. Convergence of minimizers. Now we turn our focus to the consistency of 
the regularization, i.e. the convergence of a sequence of solutions to (Pq,(k )) with 
a double-well potential approaching the double-obstacle potential in the limit as 
k —> oo, to a solution of (P<n) with the double-obstacle potential. For this purpose, 
we consider a sequence of functionals {'J'^} A , eN satisfying Assumption 3.1.2 and a 
corresponding limit functional \F 

The following theorem provides conditions under which a sequence of globally 
optimal solutions to (P^w) converge to a global solution of (P^), as k —> oo. 

Theorem 5.1 (Consistency of the regularization). Let the assumptions of The¬ 
orem 4-1 he fulfilled. The objective J : H 1 (fi) M x H 1 (fi) M x Hq a (Q; R Ar ) M_1 x 
L 2 (H;R jv ) m_1 —> R. is supposed to be upper-semicontinuous, and let N 

be a sequence of potentials satisfying Assumption 3.1.2. Assume further that T 
is given such that for every sequence {(x {k) ,y {k) )} km C H (fl) x H (H) with 

y( k ) = \]/( fc ) , (x( fc )) and (x(D,y( fc )) (a;(°°),y(°°)) strongly in H 1 (H) x H 1 (fl) it 

holds that y(°°) G 9T(x^°°^). 
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Then a sequence {(p^, ^ k \ v^ k \ u^)} keN of global solutions to (P^w) in 
H (S7) m x H (f T) M x Hq a (Q; x U a d converges to a global solution of (Pqr), 

provided that (j(<p^ k \ p^ k \ v^ k \ zt^ fc ^)} fcgN is assumed bounded, whenever U a d is un¬ 
bounded. 

Proof First note that the sequence { u ^} fc6N is bounded in the reflexive Banach 
space L 2 (fl; R' V ) M_1 . This follows either from the boundedness of the set U a d or from 
the partial coercivity of J and the boundedness of {J(p' k \ v^ k \ M^)} fceN - By 

Lemma 3.6, the sequence { (ip ^ ) yf k ) ; v ( k )) }is bounded in H 2 d jn) M xH 2 dn (n) M x 
Ho,ir(fy R Ar ) M_1 - Hence there exists a weakly convergent sequence {w/ fc! i} ;gN := 
{(<p( kl \ u( fc, ))} ZeN with limit point W:=(Jp,~p,v,u)£H^ n (Q) M x H 2 dn (Q) M x 
Hq a (Cl; R Ar ) M_1 . Moreover, since U a d is weakly closed, u belongs to U a d- 

As in the proof of Theorem 4.1, it can be shown that the limit point satisfies 
(p,J[,v) £ Sy(u). The only difference is that inclusion (4.2) follows from the above 
assumption instead of the maximal monotonicity. 

Next, we prove that ui is an optimal point of (P^). For this purpose, let (ip, //, v, u ) 
be an optimal solution of (4%-). We consider a sequence (ip ( ' k \pP z 'l) £ H^ n (Q) M x 
H 2 dn (£l) M such that 



for every <j> £ if 1 (f2) and i £ ( — 1,.., M — 2}, where v corresponds to the previously 
specified solution of (P^). Note that the operator L^ x iL 1 (H) —> H 1 (H) x 

H 1 (f1) defined by 

(5.1) L[ k \p,n) := ^-Aip+ (^o fc) ) M-M^-div(aV/j)j 


is monotone, coercive and continuous, if a £ H 2 (fl) satisfies 0 < rb\ < a(x) < 
rb 2 almost everywhere on Cl. Hence for fixed k £ N, the pair (^,- +1 ,/i) + i) of each 
subsequent time step is uniquely determined as the solution to 


(5.2) 


L 


C k ) 


(£i+i>Mi+i) = («# 


(k) 


•~(fc) Y7 —(fc) \ 

Vi - rv i+1 Vp\ ') 


where 0 < rb± < a := m((p^ kS> )h < rb 2 almost everywhere on Cl (cf. [40, Chapter 
II, Theorem 2.2]. Then, by Lemma 3.3 the sequence (<p( k > ,p,( k \v)keN is bounded 
in H 2 (Cl) M x H 2 (Cl) M x Hq a (Cl;M. N ) M ~ 1 . Consequently, there exists a subsequence 
(denoted the same) which converges weakly in the associated product space to a limit 
point (ip* , fi* ,v). In accordance with the above observations, (£>* +1 , /r* +1 ) is the unique 
solution to 

( Vi+i-Vi ^ + ( ^ + 1 v^, 0 ) - (div(m(£*)V/I* + 1 ),0> = 0, e H\n), 

(- a tp* +1 ,d>) + (d^*(p*+i)^) - (K+1,4) - (wit) = o,\/cf£ h\ci) 
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for every i £ { — 1 ,M — 2}. Note that here we also use the prerequisite that y* +1 £ 

d^oivi+i) ™ ith v (k )i +1 =y { o k) (v^i)^ Sin ^ e the 

feasibility of ((p,p,v,u) implies (c p,p,v ) £ S-^-(u), this yields (p* = (p and p* = /I 

Now we show that p^ converges strongly in (iL[) (fl)) M to p*. For this purpose, 
fix * £ {—1 ,M — 2} and define 


(5.3) 


~(fc) 
(fc) V\+1 
9i : = - 


Vl 


.(fc) 


+ 


v i+1 V(pf\ 


Vi+i - Vi 


T 


+ v i+ i\7(p*. 


By the Rellich-Kondrachov theorem g^ converges strongly in L 2 (fl) to g*. It further 
holds that g — g* = div(m(^ fc ^)V^*\) — div(m(tp?)Vp* +1 ). Hence, we have 

div(m(tp;)V(g£[ - p* +1 )) =g, {k) - g* - di.v((m($ fc) ) - m(^))V^i) =: 6^. 


Again by the Rellich-Kondrachov theorem m(<p[ k ^) converges strongly to m(tp*) in 
bF 1,5 (fI). Furthermore, is bounded in U 1 (f2). As a consequence, df ' 1 —t 0 

strongly in L 2 (fl). Applying [34, Theorem 2.3.1], we conclude 


u (k) 

Vi+l 


Vi +1 


H 2 


< c 


:(fc) 


->■ 0 . 


Next, we define £ L 2 (fl; l w ) for all i £ {0,.., M — 2} by 


qfc) PW )«i+i-PW-i)«i , j. 

u- + i := -y-hdiv(u i+ i ®p(<$_i)Ui) 


- div(u i+ i (g) P2 2 ft m(ff 1 )Vjjf ) ) 

- div(2i?(^ fc) )e(v i+ i)) - /4+iV^* 0 . 


Similarly to the proof of Theorem 4.1, it can be shown that converges strongly 
in L 2 (fl;R Ar ) M " 1 to u. 

Summarizing, the sequence {{ip( k \p( k \ v, u ^)} fegN converges towards (Ip, p, v, u) 

strongly in H 1 (fl) M x H 1 (fl) M x Hq (£7; R Ar ) M_1 x L 2 (fl; Employing the 

continuity properties of the objective functional JT, this yields 

J(Tp,p,v,u)< lim J(^ k \p^ k \v^ k \u^)< lim J(^ k \p^,v,u^) 

k—too k—too 

(5.4) < J((p,p,v,u). 


Since (<p, p, v, u ) is optimal, the assertion holds truc.D 

In summary, the optimal control problems under consideration are well-posed and 
admit globally optimal solutions. Furthermore, the chosen regularization approach is 
consistent in the sense of Theorem 5.1. 

6. Stationarity conditions. Now we turn our attention to the derivation of 
stationarity conditions for the optimal control problem. For smooth potentials ’I'o 
stationarity or first-order optimality conditions for the problem (Pip) can be derived 
by applying classical results concerning the existence of Lagrange multipliers. The 
latter approach is employed in the following theorem. 
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Theorem 6.1 (First-order optimality conditions for smooth potentials). Let J : 

xH\rt) M xiJ 0 1 CT (n ; R Ar ) M - 1 -> R be Frechet differentiable 

_2 

and let satisfy Assumption 3.1.2 such that Tq maps Hg (12) continuously Frechet- 
differentiably into L 2 (l2). Further, lefz := ((p,p,v,u) be a minimizer of (Py). Then 

_o Ai’ _ 2 Af A/T—1 

i/iere exist (p,r,q) El (12) x L (f2) x 7Jg CT (f2; R w ) , p = (p_i, ...pM- 2 ), r = 
(r_i, ...r M - 2 ), Q = (?o, ■•■9 m- 2 ), such that 


-(Pi -Pi-i) + a(m / (p i ), / Lt l+ i,p i ) - div(pi^ + i) - AVj_i 

r 

+ H'oOP*)* 7 '*-! - KU+1 - -p'(^i)Vi+l ‘ (9i+l - di) 

T 

- (p'((pi)v i+ 1 - ^y^-ra'(i ( p l )V// l+ i)(Dg l+ i) T v i+ 2 
(6.1) + 2r/(p i )e(u i+ i) : Dq l + div(p i+ i<p) 



- r*_i + 6 (m(ipi_i),pi_i) - div(^-^-^-m(pi_i)(H(?i) T ?;i + i) 
(6.2) - ft-i • V<Pi-i 

- -p(<Pj-i)(9j - 9j-i) - p(pj-i)(I?9j) T ^+i 
- (£>%_i)(p(p J _ 2 )^-i - P2 2 Pl m(<pj_ 2 )Vpj_i) 


(6.3) 

- div( 2 ?y(^_i)e(gj_i)) + p-,-_iV^_i 

(6.4) 

\ M_1 




£ [R + (27 ad - m)] + , 


/or all i = 0,..., M — 1 and j = 1, M — 1. Here, [K+(27 a[ j — u)] + denotes the polar 
cone of the set {r(w; — u)\w £ U a d Ar E M + }. Furthermore, we use the convention 
that Pi,ri,qi are equal to 0 for i > M — 1 along with g_i and ipi,fii,Vi for i > M. 
Moreover, a(f, w,p), b(fh,p), A 4 (r) £ Hg n (Cl)* are defined by (A 4 r, z) := J^fAzdx, 
(a(f,w,p),z) := J n -pdiv(fzVw)dx, (b(fh,p),z) := f Q -pdiv(mVz)dx, for func¬ 
tions f,m£ C 1 (Q),w £ £ L 2 (ffl) and z £ H 2 g n {ffl). 

Proof. Utilizing the spaces X and X and the set C given by 


X := H 2 g n ffl) M x H 2 gjn) M x x L 2 (12;R^)^ 1 , 

C :=H 2 gffn) M x H 2 gjn) M x 1^(12; x Uad, 

_2 Af _ 2 Af .. _ T A/f—'\ 

Y := {L (12))* x (L (12))* x H^(fl-R N y , 
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for tp = (tp 0 ,...,(p M -i), P = (Po,-,Pm-i), v = (v u u = (tti,u M _i) we 

define a mapping 5 : X —>• H by 

g{y,H,v,u) 

1 ( 7(^+1 ^ Pi) - div(m(<£>i)V/i i+ i) + Uj+i • V<pi ). = _ 2 
_ ( -Mi+i - + fE'o^i+i) - )^_i 

/ i(p(^)n i+ i - p(v?i_i)vi) - div( 2 r?(( / 9 i )e(w i+1 )) 

\ \ +div(n i+ i g> (p(<pi-i)vi - P2 ~ Pl pi)) - p i+1 Vp* 


M—2 


~ u >+' Ji =0 / 


Then, (P$) can be stated as min{j7"(ip, p, u, u) : (tp,p,v,u) £ C , g(tp, p,v,u) = 0}, 
with z = (<p, p 7 i>, it) an associated minimizer. The mapping <7 is continuously Frechet 
differentiable from X into Y. To see this, let us exemplarily consider the term 
div(m(<pi)V/tj+i). The other terms can be treated analogously. First note that 
div(m(<pi)V/Zi+i) equals Vm(<pi) • Vpt+i + m(tpi)Api + i where m(tpi) is given by 
m(tpi). Hence Vm(<pi) = m'(ipi)'Vipi. Assumption 2.2 implies that both superpo¬ 
sition operators ip i->- m(ip), tp i->- rn'(tp) are continuously Frechet differentiable from 
H 2 (tt) ^ L°°(n) into L°°(H) (cf. [42]). Therefore, the mappings (tp,p) —*■ m'(tp)X7tp- 
\7p : H 2 (Q) x H 2 (Q) —> L 3 (fl) and (<p,p) —> m(tp)Ap : H 2 (Q) x H 2 (Q) —> L 2 (fl), 
are continuously Frechet differentiable. This shows the continuous Frechet differen¬ 
tiability of div(m(<pi)V/Zj+i). The Frechet derivative of g in (ip,p,v,u) applied to 
(tp 5 , p 5 , v 5 , u 5 ) £ A is given by 


g'((p,p,v,u)(tp 5 ,p 5 ,v s ,u 5 ) 

f ( 7(^+1 -<Pi) - div(m'(<Pi)<p|V/i i+ i) -div(m(<pi)V/rf +1 ) \ M 2 ^ 

\ + v *+i • Vl Pi + v t +1 ' / i= _i 

( -Mi+i - A pf +1 + ^0(^+1; <Pi+i) - K( Pi )^Zl 

( \(p’{Pi)p\ v %+\ + y(p(^K 5 +i -p(Pi-i)wf) \ M_2 

+ div(u i+ i <g> (p'(<Pi-i)<Pi_ 1 v i + p(tpi-i)v 5 )) 

-div(u i+ i (g) ( P2 ~ 2 Pl m'(tpi-i)tpf_ 1 'Vpi - £2 ^ £l TO(<pj_i)V/L<f)) 

+ div(uf fl (g) (p(tpi-{)Vi - ^^-m(tpi-i)Vpi)) 

-dW(2r)'(tpi)tp 5 e(v i+1 )) - div(2 g(tpi)e(vf +1 )) 

V V -rt+iVvjf - p 5 i+1 Vtpi - uf +1 ) i=0 / 

Due to our convention for p-\ and vq , we require that p 5 _ 1 = 0 and Uq = 0. For 
the application of a result due to Zowe and Kurcyusz [46] concerning the existence 
of Lagrange multipliers, we show that g'(z) maps R+((7 — z) C X onto Y. For this 
purpose, let (0£, 0]", 0") £ Y be arbitrarily fixed. We have to show that there exists 
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a tuple (ip s , p s ,v s ,u 5 ) £ R + (C — z) such that 

^(V’f+i - Vi) ~ div(m'(p»)<^V> i+ i) - div(m(^)V/xf +1 ) 

(6-5) +u i+1 • VpJ + v 5 i+1 • Vipi = ej", 

(6.6) -/4+! - A( Pi+i - K<Pi + &o(Vi+i;Vi+i) = ©i> 

-(p'(^i)^f^+i - p'(^-i)^f-i^) + -(p(<Pi)^+i - P(<Pi-i)vf) 

T T 

+ div(u i+ i <g» {p'{Vi-i)Vi-iVi + p(<pi-i)vf)) 

-div(v i+ i ® ( P2 2 Pl m'm - P2 2 Pl m(^-i)V/zf)) 

+ div(i>f +1 (g> (p(v?»—i)i»i - P2 2 Pl m(yi_i)V/Hi)) 
-div(277'(^)<pfe(^+i)) - div( 2 ? 7 (^)e(wf +1 )) 

(6.7) - Mi + 1 V<Pi - «J+i = ©?, 

where (6.5) and ( 6 . 6 ) hold for * = — — 2 and (6.7) for all i = 0, — 1. 

As in Theorem 3.2, standard arguments show the existence of (<Po,/Zo) 6 H dn (Sl) x 

_2 

H g ^(Q) such that (6.5) and ( 6 . 6 ) are fulfilled for i = —1. Now we apply induction 
over i. Therefore, let us assume that (6.5)-(6.7) hold for i < M — 1. In order to show 
the existence of a solution to this system for * + 1 , we note that it can be written as 

^(<pf+2 - Vi+ 1) - div( to( v?i+i)V/4 +2 ) + v i +2 • v< Pi+1 = ©AO 
~Pi +2 — A( P?+2 — K( Pi+l + ^0 (Vi+2t Vi+2) = ©¥>) 
~(p(<Pi+i)vf +2 - p(<Pi)vf +1 ) + div(uf +2 <g> (p(<p*H+i - P2 2 Pl m(ipi)Vm+ 1 )) 

- div( 2 » 7 (^ i+ i)e(uf +2 )) - p 5 i+2 S/ip i+1 - u s i+2 = O v , 


for a triple ( 0 ^, 0 ^, ©«) £ (L 2 (fl))* x (L 2 (fl))* x iJg CT ( 0 ; M. N )* that only depends 
on (ip, p, v), on ip\, and vf for i < M — 1 and on (0? +1 , 0^_ 1; 0 " +1 ). But now the 
existence of a solution follows readily from Theorem 3.2 and from Lemma 3.3 when 
choosing v = p(ipi)v i+ 1 - P2 ~ Pl m(<pi)Vp l+1 as well as f 0 = p(tp i+ i), /_i = p(ipi) and 
u i +2 = 6 - Notice, here the functions p(<pi+i), m(ipi+\), rj(ipi+i) do not depend on the 
unknown <p s +2 . Further observe that we can always find a convex, affine functional 

if> : H dn (Q) R with (Di/j)z = ’Fq (p i+2 ;z) for all z £ H 2 dn (Q). Hence we deduce the 
existence of a Lagrange multiplier (p, r, q) £ Y* such that 


( 6 . 8 ) 


J'(ip,p,v,u)(ip s ,p s ,v 5 ,u s ) 


{g\Vi P, V, u)(ip 5 , p 5 , V s , u S ), (p , r, q)) 
(g'((p, p, v, u)*(p, r, q), (ip 5 , p S , v 5 ,u 5 )) 


for all (cp s ,p s ,v 5 ,u 5 ) £ H 2 g n (fl) M x H 2 dn (n) M x JTq ^.(H; R JV ) M_1 x M+(I/ ad - u). 
In order to derive the desired system for (p, r, q) from this variational equation, the 
adjoint of g'(ip, p, v, u ) has to be calculated. Exemplarily, we show this calculation for 
two terms. First, consider the term div(ui+i <E> (p' (ipi-\)ip 5 _-yVi)) which gets tested by 
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qi. Notice that for vector fields z^\ in iT 1 (0;]R Ar ) and with |gn = 0 we 
have 

(6.9) [ z {3) • div(z (2) <g> ^ (1) ) = - [ z {2) ■ {Dz {3) )z {1 \ 

Jo, Jq 

by GauB’ theorem. Hence we get 


(div(i>i+i <2> p 1 (^Pi—i)(Pi_iVi), qi) = — / v i+1 ■ (Dqi)(p'((pi~x)<Pi-iVi)dx 

J n 

= - p'■ (Dq z ) T v i+1 dx. 

J n 

Secondly, the term div(ui+i <g> — P2 ~ Pl gets tested by q j. This yields 


(div(t!j + i <g> - 


P 2 - Pi 




/ v i+ i ■ ( Dqi)( P2 Pl m((pi-i)VPi)dx 
In z 


P2 - Pi 


m(<pi- i)V/x- ■ ( Dqi ) v i+ \dx 


p S , div(- P2 Pl m(ifii-i){Dq z ) T v i+ i)dx 


since ^i-i-ilan = 0. The other terms can be treated similarly. After collecting all terms 
which contain ipf, pf and vf , respectively, it follows that 


g'((p,p,v,u)*(p, r, q) 

( ( ~Pi-i) + a(m'(tpi),pi + i,pi) - <tiv(piV i+1 ) - AVj_i \ M 1 ^ 


+'&oi l Pi)*n-i - nr i+1 - p'(tpi)v i+ 1 • i(g i+ i - i 
-{p'{(pi)v i+ 1 - £ ^ £i m'(^ i )V/r l+ i)(L)g i+ i) T v i+ 2 
\ + 2 ? 7 , ((p i )e(v i+ i) : Dqi + div(/x i+ ig l ) 

-fi-i + b(m(ipi-i),Pi-i) - dW(^^-m(ip i _ 1 )(Dq,) T v i+ i) 




iT- 0 ! 


\ 1 ■ 

/ -p(ipi-i)\(qi - qi- 1) - p(ifii-i)(Dqi) T Vi +1 

-{Dq i - 1 )(p(ip i - 2 )vi-i - £ 2^-m(ipi- 2 )\'Pi- i) 
\ - div(2?7(<pi_i)e(gi_i)) + pi-i V<pj_i 

\ ( -9,-1 ) j=1 


M—1 


Plugging this into (6.8) and using the fact that (ip 5 , p s , v 5 ,u 5 ) can be chosen arbitrarily 
in x H 2 dn {D) M x Hq j(T (H; R jv ) m_1 x M+(f7 ae i — u), we obtain the desired 

system for (p, r, q). □ 

The preceding theorem states first-order optimality conditions for problem (P^/) in 
the case of smooth double-well type potentials. In the following, we derive stationarity 
conditions for a nonsmooth potential via a limit process; compare section 7. For this 
purpose, the boundedness of the adjoint states is crucial. In order to guarantee this, 
further regularity conditions on J are required. 

Lemma 6.2. Suppose that the assumptions of Theorem 6.1 are fulfilled. Then 















M. Hintermuller, T. Keil, D. Wegner 


25 


( p,r ) £ H (12) M x H (12) M 1 and it holds that 

a(m'(tpi),iii + i,pi) = m'(ipi)Vm+i • Vft £ H 

,pi-i) = -div(m(v3i_i)Vj3i_i) £ 

-AVi-j. = -Ar<_! £ S’ 1 (12)*. 


Proof. We prove the claim by backward induction over i. For i = M — 1 we have 
Pm- l = rjf-i = 0 by convention. Now, we take the induction step from i to i — 1 
assuming that pi,ri £ H (12). This higher regularity implies for z £ 22^(12) that 


(a(m'{ipi)ifjL i+ i,Pi),z) = - pidiv(m'(tpi)zVp i+1 )dx 

■In 


= / m'{pi)zVp i+ 1 • Vpida; 

J n 


< C'||m , (^j)||L~||V/i l+ i|| L 4||Vp i || L 2||5|| i 4 

< C\\m\pi)\\ L ^\\p. l+1 \\ H 2 \\pi\\ H i\\z\\ H i 


because of V/ij+i • n = 0 on 912. Consequently, a(m'(ipi), Hi+i,pi) £ 2? (12)*. Equa¬ 
tions ( 6 . 1 ) and ( 6 . 2 ) and the assumption yield that AVj-i, , Pi-i) £ 

22 (12)*. By standard regularity arguments one shows that rj_i and Pi-i are in¬ 
deed elements of 22 *( 12 ) and the desired relations for b(m((pi-i),pi-i) and AVj_i 
follow at once. □ 

The next lemma is used in the subsequent theorem in order to prove the bound¬ 
edness of the adjoint state. 

Lemma 6.3. Let a > 0 be given and Mi and M 2 be bounded subsets of H (12)* 
and Hq ^( 12 ; M^)*, respectively. Let Ml be the set of all tuples (p, r, q;A-,h p , h r , h q \c, u; 
rh , fj, p) with 


(p, f, q) 

A 

(Jlri hp , hq ) 
(c,m) 
rh,fj,p 


£ if 1 ( 12 ) x if 1 ( 12 ) x Hi An-R N ), 

G C(H (fi); H (fi)*) be monotone, 

G Mi x Mi x M 2 , 

£ H\Ll) x H\n;R N ), 

£ L°°( 12 ) with 1/a > rh,fj > a and p>0 a.e. on 12 , 


for which the following system is satisfied: 

(6.10) — p — Ar + Ar = h r , 

T 

(6.11) —f — div(mVp) — q ■ Vc = h p , 

(6.12) -pq — div(2 rje(q)) — ( Dq)u + pVc = h q , 

T 

(6.13) - [ p\q\ 2 - {(Dq)u,q) > 0. 

t Jn 

Then the set {{p,r,q) ■ {piT,q\A\h p ,h r ,h q \c,u\m,fi, p) £ M} is bounded in H (12) x 
H\0) x 22 ^( 12 ;^). 
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In order to keep the flow of the presentation, we defer the proof to the appendix. 
Employing the preceding results, we finally perform the limit process with respect 
to the first-order optimality conditions of Theorem 6.1. 

Theorem 6.4 (Stationarity conditions). Suppose that the following assumptions 
are satisfied. 

1. J' is a bounded mapping from H (12) M x H (12) M x Hq CT (12; K Ar ) M_1 x U a d 

into the space (T? 1 )! 1 )“ xx fl; R ^)" -1 x l 2 (n ; f and 

satisfies the following weak lower-semicontinuity property 



for converging weakly in H 2 * (f 2 ) M xHo ( 12 ) M x 

x U ad to z= ( ip,p,,v,u ). 

/ \ _2 _ 

2. For every n £ N let : H d (12) —>• 1 be a convex, lower-semicontinuous 
and proper functional satisfying the assumptions of Theorem 6.1. 

3. Let «("),«(">) G x H 2 dn (n) M x Hq CT (12; R ^)" -1 x U ad 

be a minimizer for (P^w) and let {p^ n \A n \q^) € H (12) x H (12) x 
Hq 0 .( 12 ; 1^) be given as in Theorem 6.1 and Lemma 6.2. 

Then there exists an element (<p, p,v,u,p,r,q) and a subsequence denoted by 
{( 99 M, ,v^ m \ u , r, q ^)} mgN with 

weakly in weakly in Ug n ( 12 ) M_1 , 

weakly in U 2 (12;R^)"" 1 , it (m) -rn weakly in L 2 ( 12; 
p( m )^.p wea kly i n H 1 (Ll ) M , jd ™)—>r weakly in iJ 1 ( 12 ) M_1 , 
q( m ^—>q weakly in Hq. ct ( 12 ; R JV ) M ~ 1 , Tq”^ eX* weakly in H 1 ( 12 )*, 


/or all i = — 1,..., M — 2 such that for z = ( 9 ?, p, v, u) and qk := qk-i it holds that 
--(Pi ~Pi- 1) + m(9?j)'V/t i+ i ■ Pi - div(piV i+1 ) - Ar;_i 

T 

+ Af_i - nr i+ 1 - -p( 99 i)'t; i+ i ■ (<29+1 - ft) 

T 

- (p(ipi)'v i+1 - 92 Pl m'p i+1 )(Dq i+1 ) T v i+2 


(6.14) 


QJ 

+ 2rj(ipi)'e(v i+1 ) : Dq t + div(p i+1 qi) = - —( 2 ), 

Otpi 


- n -1 - div(m(^i_i)Vpi_i) - div(^—-^-m(ipi-i)(Dqi) T v i+ i) 


(6.15) 


„ dJ ( 1 


~ -p(<Pj-i)(Qj ~ Qj- 1 ) - P{<Pj-i){ D Qj) v l +1 

- {Dqj-i)(p{ipj-2)vj-i - 92 . Pl m(9? J - 2 )V/r J -i) 


(6.16) 


9J, 


-div(2»y(</j j _ 1 )e(g J -_i))+Pj-_iV^_i = ^-(z), 

dJ, 


-^■{z)-qG [R+(f/ Q d-u)] 


(6.17) 
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Proof. 1. In the first step, we show the boundedness of r^ n \ q^)} N 

_ 1 M _ 1 M A/i _-I 

in H (fL) x H (fl) x Hg a (fl;S. N ) . Moreover, the boundedness of the se¬ 


quence {(¥> (n) ,/r (n) ,u( n ),u( n ))} ngN in H dn (fl) M x H dn (fl) M x H^fl-, R JV ) M " 1 x 
L 2 (fl; R jv ) m_1 follows from Lemma 3.6. For i = 0, ...,M— 1, j = 1 1 and 

n £ N the adjoint system for (P^(nj) corresponding to (6.1)-(6.4) can be rewritten as 


(6.18) 


^S-Arn+^r (vrrr-n = 

“4-1 - div(m(^” ) 1 )V^” ) 1 ) - = B©^, 

- div(2r ? (^"\)e(g}"\)) +pj” ) 1 V^” ) 1 

(Dqj—l')(p(Pj 1 —2 )^j'—~ ^ 2 Pl ^(^- 2 )V^" } i) = ©,"•_!, 
where the functionals b£"\ Bp"" 1 and Bg""* are given by 


J") 


J") _ o(™) 


0 


(ra) 


Oipi T 


™'( 4 n) ) v 4 +i -p[ n) - div(pf ) ^" ) 1 ) - Kr^ 


- ^P(^ n) ) 4+1 • (?l+i - 4" } ) + 2 »?(^i n) ) e(4+i) : + div(/i J ( ” ) 1 ^ n) ) 


- o^n - 

dj, 


2 

P 2 - Pi 


■ m '{Pi )V/x^"i)(-Dg'+i) T uj" 2 

TO (^-i)(- D 9i n) ) T 4+i). 


e ^-'= aS ( " l ” >) + div( - 
ejt, = g( z <">) + ^(di)4 ( " > -/-(d- > .)(04 < " ) ) T *'S- 

Here, z ©- 1 denotes the tuple (^ n \ //("),«("), u(" i). We prove the boundedness of 
{(p("),r("),g , ("))} ngN in iJ 1 (f2) x H 1 (fl) x Hg a (fl;M. N ) M 1 by backward induc¬ 
tion over i. If i > M — 1, then (p[ n \ r[ n \ q^) = 0 by convention. In the induction 
step assume that for i £ {0,..., M — 1} and for j >f the sequence j (p^, r < 7 ^ ) j 

is bounded in H (fl) x H 1 (fl) x Hg a (fl; 1^). This and the assumption on J imply 
that {(Bp"^!, Brli-nB^-i)} is bounded in (T? (fl) x Tf\fl) x Hg cr (fl- 1 R N ))*. 

To see this, we exemplarily consider first 2r]((p\ n ^) e(v^f):Dq^ n \ which is bounded by 
ll 2 ^ 70 ) e(4+i) : D Qi n) \\LVB < C\\p(p[ n) ) 11l=o 1| e(4" 1 )11ls11 Dq\ n )11p, 2 

< C\\V(p[ n) ) 

and secondly — P 2 ~ pl m!(ip^) V (Dq^ ± ) T v- ” 2 , which we bounded using 

■ I P2 Pl /, (n)\r-> (n) / n (n) \T (n) i 
I-o- 171 Wi Pi+l( D( li+l) v \+?\\l*/s 


< C||- 


P 2 


ra' ( < P l (n) ) 11 L“ 1IVp ■ + ! 11 L» 11 £>4+1 11 lA 14+2 11 + 6 


< C'll - P2 Pl m / (4” ) )lk<”|l4+il|gHl4+illg 1 ll4+2llg 2 - 


2 
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Consequently, these terms define continuous linear functionals on H (fl), that are 
bounded independently of n. The other summands can be estimated similarly. 

In case of i > 0 we apply Lemma 6.3 to 


(p, r, q; A; h p , h r , h q \c, u; m , r), p) 

d-n/Kd-l)*!-! - — - — 7 .'' C-'i i i j’V-'i 


Note that due to divt/"^ 


= 0 we have 


div u = p{4-\) v4i • - div( 


(n) 


, p2 - Pi 


m{4- 2 )^Pi-i) 


(n) 


P 2 - Pi 


4 -\ ' v< P;-i - div(m(^” ; 2 )V^” j 1 ) 


,(«) 


(n) 




1 P2 Pi ^ (n) (n) ^ ( (n) ^ t (ra) ^ 

-~-(<Pi-l - Vi-2) =- {Pm- 1) - P{Vi—2 ) )' 

T Z T 


With the help of J n ((Dq)u, q) = — f Q q • div(g ® u) (cf. (6.9)), (3.7) yields 


1 

r 


p\q\ 2 dx - (( Dq)u,q) 


- [ p(4-i)\qt\\ 2 - 1(P(4- 1 ) - p{4^2))\44i\ 2 dx 
T Jn 2 

[ {p{4-\) + p{4-i))\4-\\ 2 dx > 0 , 

J Cl 


because of p(4- 2 ) — 0 an d p{4- 1 ) — 0 almost everywhere. Hence Lemma 6.3 implies 

the boundedness of (p-"\,r^, q4i) in -ff^fi) x iL 1 (H) x 

The case i = 0 needs some modifications in order to be treated by Lemma 6.3 
since (6.3) is not defined for i = 0. In this case we set (g, h q , c, u, f))={ 0, 0, 0, 0, p{44)) 
together with the definition of the remaining quantities as in the case * > 0. Now, by 
Lemma 6.3 we conclude the boundedness of {44i , 44 ) in H (fi) x H (H). Moreover, 

from (6.18) it follows that also (’Lp"'* (44)* r 4 i) remains bounded in if^fi)*. 

2. With the bounds derived in step 1 and with the usual compact embeddings of 
Sobolev spaces, we can pass to a subsequence with the desired convergence properties. 

3. Now we pass to the limit in the the adjoint systems corresponding to (6.1)-(6.4) 
for (P\j/(n)). The limits for the equations (6.1) and (6.2) are considered in H (H)* 
and the limit for (6.3) in #q ^(fl; R^)*. In the linear terms we can pass to the 

limit at once. For m!{44 /4+i '44 we have that m'(44 ) converges strongly in 
L°°(fi) to V pd+i strongly in L 6 ~ s (fl) to V/q+i and p- n) weakly in L e (fl) 

to pi . Hence, m' (44)^ 4i +1 ' p4 converges weakly in H (H)* to m'(tfi)X7p,i + 1 • pi. 

For P2 ~ Pl m!{44 ) Vp 4i (-Pgi+i ) T ^~+2 we note that P2 ~ Pl m'(44) anc l 4+2 converge 

strongly in L°°(H) to P2 ~ pl m'(pi) respectively Vi+ 2 , V/4+i strongly in L 6_e (H) to 
V/q +1 and Dq4 weakly in L 2 (fl) to Dqi. Therefore pa ^ Pl ml{44 ) V/4 +1 (-P<?j+i 44+2 
converges weakly in IL 1 (H)* to P2 ~ Pl m' (tpi)V p,i+i(Dqi + i) T rg +2 . 

For div( P2 ~ Pl m(44i){^44 4 v 4 \) we use that P2 ~ Pl m(44i) an d v 4i converge 

strongly in L°°(H) to P2 2 Pl m(ipi- 1 ) respectively Vi+ 1 , and q4 weakly in L 6 (H) to 
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qi. As a consequence, div( P2 ~ Pl m(<p ) (Dq^ ) T t>|”\) converges weakly in H 1 (LI)* 
to div( P2 ~ Pl m(yj,;_i)(-Dgi) T t)i + i). For the convergence of div(2r](ip[ r f > 1 )e(qlZ > 1 )) note 
that r) (</j-"\) converges strongly in L°°(Ll) to and weakly in P 2 (fl) to 

e(qi-i). Hence, div(2r?(^"\)e(g-"\)) converges weakly in Hq a (Lt; R w )* to the limit 

div(2?7((^i_i)e(gi_i)). Apart from 'Fq"' ) (<^™' ) )*t^ 1 , all remaining terms appearing 
on the left hand sides can be treated similarly. Moreover, our assumptions on J 
imply that J'(^ n \ H^ n \ v^ n \ u^) converges weakly to J' (ip,\i,v,u) in (H 2 g n (Ll) M x 
H 2 dn (n) M x x L^R^" 1 )*. 

Consequently, by (6.1) also converges weakly in to 

some Aj_i. Therefore, we arrive at the system (6.14)-(6.16). Finally, notice that 
for all y £ U a d and with := {^ n \ ^ n \v^ n \u^) and := q^ 1 by the weak 
lower-senricontinuity of ^ an d the weak and strong convergence of the sequences 
involved we deduce that 

(|iSr(z) - Q,y~u ) = (§£ (z), y) - (§£(2), u) — (q,y — u) 

> liminf (<^(^")),y) - (|( Z W),«(")) - (q {n \y- « (n >)) 

= liminf (^(^ (n) ) - q {n \y~ w (n) }) 

> 0 , 

due to the optimality of z^ for (P^(„)). This shows (6.17) and finishes the proof. □ 
Remark 6.1. We point out that a tracking-type functional, like, e.g., 

1 2 C 2 

J((p,p,v,u) := - \\ipM-i ~ <Pd\\ + 2 I|m||(l2 )(M -i) , £ > 0 , 

with (fid £ L 2 (Ll) a desired final state, satisfies the assumptions of Theorem 6.4-. 

Remark 6 .2. If the set U a d is bounded, Theorem 6.4 holds also true for a sequence 
{()))neN °f stationary points for (P^m )■ If it is unbounded, then 
the result can still be transferred to sequences of stationary points by assuming that 
the sequence { M ^} ngN is bounded in L 2 ( f2;R Ar ) M_1 . 

7. Stationarity conditions in case of the double-obstacle potential. In 

this section, we apply the developed theory to the initially stated optimal control 
problem associated to the double-obstacle potential. For this purpose, let ifo be 
defined as in Assumption 3.1.1 and set 7 := dif 0 C R x R. Then we define the 
sequence of approximating double-well type potentials as follows. 

Definition 7.1. Let a mollifier <4 £ C^R) with supp£ C [—1,1], = 1 and 

0 < £ < 1 a.e. on R, and a function 6 : R + —> R + , with 0(a) > 0 and -> 0 as 
a —>• 0, be given. For the Yosida approximation 7 a with parameter a > 0 of 7 define 

C a(s) ■■= -((-), 7a := 7a *Cfl(a), 4>0 a ( s ) : = 7a(0 dt, 

a \a/ J 0 

'J'oa(c) := / C^Oa ° c)(t)dt. 

■hi 

Moreover, we set a n := n _1 , 4 'q 1 '* := fo an . 
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Remark 7.1. We note that '&q J ^ can be identified with the superposition oper- 

_ _ _ 2 

ator corresponding to 7 a „, cf. [25]. Since 7 ' is bounded and since Ftg n (Pl) embeds 

continuously into L 2 ~ s (Pl) for S > 0, it follows that maps Hg (PI) continuously 

Frechet-differentiably into L 2 (Q), see, e.g., [20]. 

In order to obtain a stationarity condition for the optimal control problem of 
CHNS with the double-obstacle potential we pass to the limit (with the Yosida param¬ 
eter) in a sequence of optimal control problems with approximating double-wcll-type 
potentials. 

Theorem 7.2 (Limiting e-almost C-stationarity). Let 'I','”' 1 . n £ N be the func¬ 
tionals of Definition 7.1, and let the tuples p^ m \v^ m \u^ m \p^ m \r^ m \q^ m ^), 

(ip, p,v,u,p,r,q) and J be as in Theorem 6 . 4 . Moreover, let A : R —> K be a Lips- 
chitz function with A(ipi) = A(^ 2 ) = 0. For 


a „■ 


(m) 


:= 




A<“> 


:=T, 






for i = 0,..., M, and for cq denoting the limit of a[ m \ it holds that 

{ai,A(<fi) ) L 2 = 0 , = 0 , 

{ai,n -1 ) L 2 = 0, liminf( ) L 2 > 0. 


Moreover, for every e > 0 there exist a measurable subset Mf of Mi := {x £ PI : ip 1 < 
<Pi(x) < 1 P 2 } with \Mi \ Mf | < e and 

(Xi,v) =0 Vu £ Ft (12), u|o\Mf = 0. 


Proof. 1. The subdifferential 7 satisfies yA(x) = 0 if (x, y) £ 7 . Since (<£*, a*) £7 
a.e. on 12 and since a* £ L 2 (Pl), integration yields the complementarity condition 

( ai,A(ipi) ) L 2 = 0 . 

2. Now we show that ( Xi, A(ipf) ) L 2 = 0. It is well-known that the superposition 
Pk of the metric projection px of R onto K := [ipi, 1 P 2 ] maps Ft (12) continuously into 
itself. Denoting by La the Lipschitz constant of A, it holds that |A(s)| < LAmin(|s — 
ipi\, |s — ip 2 \) for s £ R. Using |%(s)| < A for all s and j' a ( s ) = 0 for ipi + 0(a) < s < 
ip 2 — 0(a) (cf. [25]) yields 

\(X ( r\A(P K (p[ m) )))L^\ 2 = |(r^,v I ;W"(^) )A( p if( ^) )))i2 |2 

< Wr^Wh [ |7' m (^ m) )A(Px(^ m) ))| 2 

Jn 

. f 1^1 |i (m) i| T @((%m)\ n 

< ( 1121 11 r) ’\\ L 2 La- -- -»• 0 

\ Q-m J 

as m —> 00 and consequently 

= lim( Al” 1 , A(f K (Y>‘ m, » ) l . + lim(A<”>, A( v <”>) - A(iWrf">)))*. (n) 

= 0 , 

which implies (Xi, A(ipi)) = 0 since (p[ m ' > converges strongly to ipi = Px(fi) in H 1 (Pi). 
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3. Denoting g m (s) := 7 Qm (s) - 7a m ( s ) 7r ( s ) with s ~Pk(s) —■ tt(s) yields 




(r, ( "J. „ + (A'”>, - PkW?" )) v . 


Since |g m (s)| = | 7 « m (s) — 7a m (s)7r(s)| < C e ^ m ' ) for m sufficiently large (cf. Lemma 
4.2 in [25]), the first term on the right-hand side converges to 0 and the second one 
as well because of the strong convergence of and (Pk (^)) to tpi in iL 1 (D), 

respectively. 

4. The property liminf(A )^2 > 0 follows readily from the monotonicity 


of 


(Pi )■ 


5 . The convergence properties of imply that the subset G := {x G fl : 
—> (fii(x) asm 4 00 } of Q has full measure (i.e. |G| = |D|). Therefore, for 
every ifGnMjwe can find m 0 (x) G N with ipi + 9(a m ) < ip[ m \x) < ip 2 — 9(a m ) 
for all m > mo{x). Thus, (x) = 7 ^ m ( x))r[ m ' > ( x) converges to 0 on G PI Mi. 

Using Egorov’s theorem shows that for every e > 0 there exists a subset Mf of GflAf; 
with |Mi \ Mf | < e such that \[ m ' > converges uniformly to zero on Mf. Hence, we 

obtain (A i,v) = lim(A- m \w) = 0 for every v G H l (D) with w|n\M? = 0. □ 

In combination with the results from Theorem 6.4, Theorem 7.2 states stationarity 
conditions corresponding to a function space version of C-stationarity for MPECs, cf. 
[22, 23], 


8 . Conclusion. Our specific semi-discretization in time for the coupled CHNS 
system with non-matched fluid densities represents a first step towards a numerical 
investigation/realization of the problem. Most importantly, it preserves the strong 
coupling of the Cahn-Hilliard and Navier-Stokes system which, in the case of non- 
matched densities, is additionally enforced through the presence of the relative flux 
J. As a result, wcll-posedness of the time discrete scheme is guaranteed and energy 
estimates mirroring the physical fact of decreasing energies can be argued. Such an 
energy property is not clear for the time continuous problem at this point in time and 
might be the subject of further research. 

Concerning the potential chosen in the Ginzburg-Landau energy, we note that 
while the existence of global solutions to the optimal control problem can be shown 
for both cases (i.e., for double-well and double obstacle potentials) simultaneously, 
the derivation of stationarity conditions is more delicate. In fact, the double-obstacle 
potential gives rise to a degenerate constraint system with the overall problem falling 
into the realm of mathematical programs with equilibrium constraints (MPECs). In 
our approach, the constraint degeneracy is handled by a Moreau-Yosida regulariza¬ 
tion approach (resulting in an approximating sequence of double-well-type potentials) 
and a subsequent limiting process leading to a function space version of so-called C- 
stationarity. For the underlying problem class, our limiting version of C-stationarity 
is currently the most (and, to the best of our knowledge, only) selective stationarity 
system available. As an alternative analytical approach, one may want to pursue set¬ 
valued analysis in order to derive stationarity conditions directly, i.e., from applying 
variational geometry (contingent, critical and normal cones) and generalized differen¬ 
tiation. This, however, is usually not possible by simple application of available tools, 
but rather by expanding current technology. It, thus, may serve as a subject of our 
future work on this problem class. 
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Finally, we point out that the constructive nature of our derivation of stationar- 
ity conditions facilitates a numerical implementation of the approach which can be 
exploited in future investigations of these problem types, both, from a numerical, as 
well as, from a practical point of view. In [26], this has already been effectively done 
for the case of matched densities. 

Appendix A. Proof of Lemma 3.4. 

Proof. Let L e : H 1 ( f2) —> H 1 (fl) be defined by 

(L e {p), 4>) ■■= (-Ap,< p) 


(A.l) 


- (g x + max(— 0 i, O)0 £ (<p - 

-ifx) 

+ min(- 

-51,0)06(^5 

where <f> 

£ h\ci) 

and 0 £ is defined by 






f 1 

if x 

< 0 , 




0 s (x):={ 1 -f 

if 0 

< X < £. 

i 



l 0 

if x 

> £. 


Since g\ 

£ L 2 (n) 

and 9 e (<p - if i), 9 e (if 2 - p] 

)£L' 

“(«), it 

holds that 


(A.2) \\ 9 i+^x(-g 1 , 0 )e e ((p-ifi) + mm(-g 1 ,0)e s (ip- ^ 2 )|| < llffill ■ 

We show that for every 0 < e < min(— ipi,ip 2 ) there exists a unique tp e £ n K 
such that 

(A.3) L e {y e ) = 0, 

In fact, for every w,v € H 1 ( fl), it can be seen that 

{L e (w) — L s (v), w — v) > / |Vw — Wv \ 2 dx 

J n 

where we use the monotonicity of d E . By Poincare’s inequality there exists a constant 
C > 0 such that 

(L e (w) - L e (v), w - v) > || Vw - Veil 2 >C\\w- vf H1 . 

Consequently, L e is strongly monotone and coercive. Since L e is also continuous 
on finite dimensional subspaces of H (f2), [31, III: Corollary 1.8] is applicable which 
yields the existence of ip e £ H (fi) with L e (ip e ) = 0. 

Due to the definition of L e and inequality (A.2), we have A ip e £ L 2 (Q). By [34, 
Theorem 2.3.6] and [34, Theorem 2.3.1] there exists a constant C i > 0 such that 

(A.4) \\Ve\\ H * < Cl ||A<p e || + ||<p e ||. 

In combination with (A.2) and Poincare’s inequality, this leads to 

(A.5) \We\\ H 2 < C 2 ||si||. 

Now, we set /3 £ := ip e — min(<p E , iff) > 0 and observe that 

l|V/3 £ || 2 = [ 

Jn 


(A.6) 


V(vj e - i/j 2 )Vf3 e dx = {—Aip e , f3 e ) 


M. Hintermuller, T. Keil, D. Wegner 


33 


where fix := {a; S D : /3 E (x) >0} = {x £ 17 : <p E (x) > fa>ip l + £ }- By equation 
(A.l) and (A.3), this leads to 


liv/y 2 


(fli + max(- 5 i,O) 0 e (<£ e 


fa) + min(- 5 i, 0 ) 0 e (fa 



+ min(— gi, 0))/3 e dx < 0. 


Ve))Pedx 


Thus, /3 e = 0 and therefore tp E < fa almost everywhere in f2. 

In a similar way, we prove that tp e — max.(ip e ,fa) = 0 and therefore tp E > fa 

— _ 2 

almost everywhere on fl. Hence <p E is contained in H (f l) flK. By inequality (A.5), 

_ 2 

the sequence {<p e }e_s.o is bounded in H (fl) and there exists a weakly convergent 
subsequence (denoted the same) such that p e —^2 fa with ||(p *|| ff2 < C 2 ||< 71 1|. Since 
K is weakly closed, it contains fa. 

For arbitrarily small 0 < <5 < min(— fa, fa), let v £ K be such that fa + 5 < v < 
fa — 5 almost everywhere in fl. Using equation (A.3) and the monotonicity of L e , we 
infer 

0 < {L e (v),v - ip e ) = (-Av,v - (fie) - / (gi+max.(-g 1 ,0)9 E (v-fa) 

Jn 

+ min(—gi, 0)0 E (fa - v))(v - tp E )dx 

= (-A v,v-<p E )- / gx(v-ip E )dx 
Jn 

for every 0 < e < 5. For e —> 0 this leads to 


0 < (— Av, v — fa) — I g\{y — fa)dx. 

Jn 

Since 6 > 0 can be chosen arbitrarily small, the last relation holds for every v £ IK 
via a limiting process. Applying [31, III: Lemma 1.5] once more, this implies 


0 < (—Afa,v — fa) — f gi(v — fa)dx 1 Vu £ K. 

Jn 

Due to the uniqueness of the solution for our variational inequality problem, this 
yields the assertion.^ 

Appendix B. Proof of Lemma 6.3. 

Proof. Testing (6.10)-(6.12) by tt , p and q, respectively, and summing up we get 

T{h r ,f) + (h p ,p) + (h q ,q) 

= r(Vf, Vr) + r(Af, f) + (mVp, Vp) 

+ ~(pq,q) ~ (( Dq)u,q) + ( 2 fje(q),e(q)) 

T 

- T ll f 11 ^( 0 ) + ^(n) + ll«lllri„(n^)) 

for a positive constant C depending only on a and on the constants in Korn’s and 
Poincare’s inequalities. This estimate yields the assertion. □ 
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